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ABSTRACT 

We measure the redshift-dependent luminosity function and the comoving radial den- 
sity of galaxies in the Sloan Digital Sky Survey Data Release 1 (SDSS DR1). Both 
measurements indicate that the apparent number density of bright galaxies increases 
by a factor w 3 as redshift increases from z = to z — 0.3. This result is robust 
to the assumed cosmology, to the details of the if-correction and to direction on the 
sky. These observations are most naturally explained by significant evolution in the 
luminosity and/or number density of galaxies at redshifts z < 0.3. Such evolution 
is also consistent with the steep number-magnitude counts seen in the APM Galaxy 
Survey, without the need to invoke a local underdensity in the galaxy distribution or 
magnitude scale errors. 
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1 INTRODUCTION 

Measurements of the galaxy luminosity function (LF) and 
its evolution provide important constraints on theories of 
galaxy formation and evolution. It is currently believed 
that galaxies formed hierarchically from the merger of sub- 
clumps, with the peak of star format i on ra te occurring 
around redshifts z ta 2-4, ee. lCole et alJ i2000Tl . Since then, 
galaxies are thought to have evolved mostly passively as 
their stellar populations age, with occasional activity trig- 
gered by interactions with other galaxies. 

Significant evolutio n in the LF ha s been measured 
since redshift z ~ 1 fee. iLillv et"aflll995t lEllis et al.lll996t 
I Wolf et al J 120051) . but most existing galaxy samples have 
been too small to directly constrain evolution at more re- 
cent epochs. By combining three different redshift surveys, 
lEalesl <ll993l) was able to demonstrate that the amplitude 
of the LF increases by a factor » 3 in the redshift range 
< z < 0.4. He also found evidence for an increase in the 
amplitude of the LF to z « 0.15-0.2, although the strength 
of this low-redshift evolution was poorly constrained. 

The Sloan Digital Sky Survey (SDSS. lYork et al.l200Cft . 
provides an ideal sample with which to measure galaxy evo- 
lution at low red shifts. A recent determ ination of the LF at 
redshift z = 0. 1 jBlanton et alj|2003d) . allowing for galaxy 
evolution, showed that the typical r-band galaxy luminosity 
brightens by ~ 1.6 mag per unit redshift. Here we further 
investigate evolution in the r-band LF of SDSS galaxies. We 
describe the data sample in Section and estimate the LF 
in four redshift slices in Sectional The effects of galaxy evo- 
lution are investigated in an alternative way in Section [I] 
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where we estimate the radial density of galaxies. Compar- 
isons with observed number-magnitude counts are made in 
Section |5] and we conclude in Section 



2 DATA SAMPLE 

We use galaxies from the Sloan Digita l Sky Survey Data 
Release 1 (DR1, lAbazaiian et alJl2003l) . The SDSS is per- 
forming five ban d CCD imagin g over an area ~ 1 0, 000 
square degrees jFukugitaetal 1 Il996t lc unnetail Il998 ~ 

M 



iHogg et alJ200ll . lSmith et aJ2002,|Pier et al .120031) . A mul- 
tifibre spectrograph is measuring spectra and redshifts for 
a subset of sources detected in the imaging data; here 
we consider only ga laxies in the main flux-limited sample 
JStrauss et al J 12002ft. A techn ical summary of the survey 
is given m lYork et alJ J2000l) and a description of SDSS 
photometric and spect roscopic parameters may be found in 
IStouehton et alJ feOOSft . 

The magnitude limit of the main galaxy 
survey has been set at an extinction-corrected 
iSchlegel. Finkbeiner. fc Davisl Il998l) r-band Petrosian 
magnitude r < 17.77. This magnitude limit was chosen as 
test year data demonstrated that it corresponds closely 
to the de sired target density of 90 objects per square 
degree, see lStrauss et al] i2002ft for full details of the target 
selection algorithm. However, a small amount of early data 
included in DR1 was taken when the magnitude limit was 
set at r < 17.6, and so we adopt this brighter magnitude 
limit here. We thus select galaxies from the DR1 catalogue 
with extinction-corrected r-band Petrosian magnitude 
r < 17.6 and with one or more of the TARGET_GALAXY, 
TARGET_GALAXY_BIG or TARGET_GALAXY_BRIGHT_CORE bits 
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set in the primTarget bit mask; see lStouehton et al.1 (120021) 
for a description of these target classifications. 

These selection criteria yield a sample of 162,989 target 
galaxies over 2099 square degrees. Of these target galaxies, 
91,921 have had a spectrum observed and 91,611 have a 
redshift measured with a confidence of 80% or higher. The 
sampling rate for our sample, defined as the number of galax- 
ies with reliable redshifts divided by the number of target 
galaxies, is thus / = 0.562. This low sampling rate is sim- 
ply due to the fact that SDSS spectroscopic observations lag 
imaging observations; the effective spectroscopic area is 1360 
square degrees. The fraction of targets observed spectroscop- 
ically is further redu ced by limitations on the placement of 
spectroscopic fibres dBlanton et al.ll2003af) . The fraction of 
spectra yielding reliable redshifts is more than 99.6%, with 
no discernible dependence on apparent magnitude. Whilst 
we naturally cannot demonstrate that spectroscopic success 
rate is independent of redshift, the fraction of galaxies with- 
out reliable redshifts is only 0.4%, and so spectroscopic in- 
completeness will have a negligible effect on estimates of 
galaxy evolution. 

Heliocentric velocities are converted to the Galactocen- 
tric frame using i>Gai = VHelio + 220 sin(i) cos(&). Unless oth- 
erwise stated, we assume a Hubble constant of Ho = 100 
km/s/Mpc and an Qm = 0.3, £Ia = 0.7 cosmology in calcu- 
lating distances, comoving volumes and luminosities. 

When estimating intrinsic galaxy luminosities, it is nec- 
essary to correct for the fact that a fixed observed passband 
corresponds to a different range of wavelengths in the rest 
frames of galaxies at different redshifts, the so-called K- 
correction. The it-correction depends on the passband used, 
the redshift of the galaxy and its spect ral energy distribution 
(SED). Here we use kcorrect v3_lb felanton et al.ll2003bl) 
in order to estimate and apply these corrections. Briefly, this 
code estimates an SED for each galaxy by finding the non- 
negative, linear combination of three template spectra that 
gives the best-fit to the five SDSS magnitudes of that galaxy. 
Rather than estimating luminosities in the rest-frame of 
each galaxy, we use kcorrect to estimate l uminosities in 
a passb and blue-shifted by z — 0.1. Following iBlanton et alJ 
(2003b) we denote the r-band in this frame as 01 r. The 
advantage of this choice of restframe is that galaxies at red- 
shift z — 0.1 (close to the mean for the SDSS main galaxy 
sample) have /^-corrections independent of galaxy type, and 
these corrections are on average smaller in amplitude than 
^-corrections at redshift zero. 



3.1 Shape 

The probability of seeing a galaxy of luminosity Li at red- 
shift Zi in a flux-limited catalogue is given by 



Pi cc cf)(Li 
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cj>(L)dL 



(1) 



where L m in(z;) and L ula ,^(zi) are the minimum and maxi- 
mum luminosities observable at redshift Zi in a flux limited 
sample. In the STY method, the likelihood C = Y[Pi (where 
the product extends over all galaxies in the sample) is max- 
imized with respect to a set of parameters describing the 
function 4>( L). For example, if we assume that <f>(L) is de- 
scribed by a lSchechterl dl976l) function, 

0(L)dL = 0*(il) Q exp(-il)d(i 7 ) 5 ( 2 ) 

we maximize the likelihood with respect to a and L* . Errors 
in the Schechter parameters are estimated by the jackknife 
method, whereby we subdivide the galaxies into 20 roughly 
equal sub-samples and estimate the Schechter parameters 
omitting each sub-sample in turn. The variance in parameter 
x is then given by 
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where N = 20 is the number of sub-samples and x is the 
mean of x. 

For the highest redshift sub-sample considered here, 
z > 0.2, only galaxies more luminous than Mo.i r ~ —21.5 
make it into the SDSS main galaxy sample. We thus have 
only a range of about two magnitudes over which to fit the 
Schechter function. Consequently, the accuracy with which 
the Schechter parameters can be determined is limited, par- 
ticularly the value of the faint-end slope a. We therefore also 
estimate the luminosity function using the SWML method 
of EEP in which <f>(L) is parameterized as a set of numbers 
(f>k in equally spaced magnitude bins. The likelihood C is 
maximized with respect to 4>k applying constraints as de- 
scribed in EEP. Also following EEP, errors on the (f>k are 
estimated from the information matrix. 



3.2 Normalization 

We use the following estimator of the space density n of 
galaxies: 



3 GALAXY LUMINOSITY FUNCTION 



We es timate the LF using the lSandaee. Tammann fc Yahil| 
lll979l STY) parametric maximum likelihood method and 
the stepwise maximum likelihoo d (SWML) method of 
lEfstathiou. Ellis fc Peterson! jl98St EEP). These estimators 
are unbiased by density inhomogeneities and have well- 
defined error properties. Both methods assume that 4>(L) 
has a universal form, i.e. the number density of galaxies is 
separable into a function of luminosity times a function of 
position: n(L, x) = 0(L)p(x). Using these estimators, the 
shape of <fi(L) is determined independently of its normaliza- 
tion. 



fn = ^ w(zi) 



dVS(z)w(z) , 



(4) 



where / is the sampling rate, S(z) the galaxy selection func- 
tion and w(z) a weighting function. The selection function 
for galaxies with luminosities L\ to L2 is 

S(z) = / cf>{L)dL / / cf>{L)dL . (5) 

-/ max(L m j n (2),ii) / J h\ 

Note that the integration limits in the numerator depend on 
the A'-correction. In this case, we estimate the if-correction 
using an SED created from the mean template coefficients 
of all galaxies in the sample. 

We adopt the weighting function 
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Figure 1. Distribution of the V/V ma x statistic for all main sam- 
ple galaxies in SDSS DM. 

1 f r ° 
w{z) = J 3 (r c ) = / r 2 £(r)dr (6) 

[1 + 47r/nJ 3 (r c )S'(2)] _/„ 

where £(r) is the two point galaxy correlation function. Pro- 
vided J 3 (r c ) converges on a scale r c much smaller than the 
depth of the survey, then the weighting sch eme (eq.El min- 
imize s the variance in the estimate of n llDavis fc Huchral 
1982). Larger values of J 3 weight galaxies at high red- 
shift more highly; we adopt 4ttJ 3 » 32, 000/i~ 3 Mpc 3 . This 
value comes fro m integrating the tw o-point galaxy correla- 
tion function of IZehavi et~al1 (I2002I) . £(r) = (r/6.14)" 175 , 
to r c = 50/i~ Mpc; at larger separations the value of J 3 be- 
comes uncertain. However, the results are not too sensitive 
to the value of J 3 , the estimated density decreasing by 7% 
if J 3 is halved. 

When normalizing the non-parametric SWML estimate, 
the integrals in © are replaced by sums over magnitude 
bins, with appropriate weighting of partial bins in the nu- 
merator. 

3.3 Results 

3.3.1 Full sample 

We first estimate the luminosity function for all 90,275 
galaxies in the SDSS DR1 with redshifts greater than 0.001 
and with absolute magnitudes in the 01 r band between 
—24 and —16. The first indication of significant evolution 
in this sample comes from the skewed distribution of the 
ISchmidtl (fl968 : 1 V/Vmax statistic (Figure[TJl. We find a mean 
{V/Vmux) = 0.523, whereas one would expect (V/Vmax) = 
0.500±0.001 for a homogeneous distribution of 90,275 galax- 
ies. 

The luminosity function of this full sample is shown in 
Figure|21 A Schechter function with parameters a = — 1.23± 
0.02, M* = -20.63±0.02 and (f>* = 0.0194±0.0010/i 3 Mpc -3 
provides good agreement with the non-parametric SWML 
estimator except for the bright end, M r < —22.5, where 
we see in the SWML estimate a higher density of galaxies 
than predicted by the Schechter function. This bright-end 
excess can be explained by evolution of the LF: the most lu- 
minous (rare) galaxies are likely to be seen at high redshift 
and, as we shall see below, high-redshift galaxies have an en- 
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Figure 2. Luminosity function for all main sample galaxies in 
SDSS DM. Symbols denote SWML estimates (the error bars are 
smaller than the symbol size for most of the points) and the curve 
shows the best-fit Schechter function. Note that no allowance for 
evolution has been made in this estimate. 



Table 1. Number density n of galaxies in the range —24 < 
Mo.i r < —21.5, in units of 10~ 4 /i 3 Mpc — 3 . For comparison, bb 
shows galaxy density inferred from the evolving Schechter LF of 
iBlanton etail 12003 d) . 



z 


^gal 


n 


TLB 


0.078 


750 


2.48 ±0.16 


2.54 ±0.07 


0.130 


2477 


3.64 ±0.15 


3.33 ±0.09 


0.178 


5559 


4.60 ±0.14 


4.21 ±0.11 


0.226 


3325 


7.36 ±0.19 


5.26 ±0.14 



hanced luminosity or density relative to low-redshift galax- 
ies. The Schechter function fit will be mostly constrained by 
galaxies at an intermediate redshift (z ~ 0.1). Additionally, 
it is possible that some galaxies with extreme luminosities, 
Mo.i r < — 23, may have incorrectly measured fluxes — see 
below. 

Note that any evolution in the luminosity function will 
render invalid our assumption that the number density of 
galaxies can be separated into a function of luminosity 
times a function of position. In order to obtain an unbi- 
ased estimate of the LF, one needs to allow for evolution. 
IBlanton et ail {2003c) have already shown that the r-band 
luminosity function of galaxies in the SDSS can be described 
by a Schechter function whose characteristic luminosity L* 
brightens by ~ 1.6 magnitudes per unit redshift. Here we 
investigate the evolution of the galaxy LF in an alternative 
way, simply by subdividing the sample into slices in redshift. 

3.3.2 Redshift slices 

Figure|3]plots the r-band luminosity function for galaxies se- 
lected in four redshift slices: 0.001 < z < 0.1, 0.1 < z < 0.15, 
0.15 < z < 0.2 and 0.2 < z < 0.3. The points with error 
bars show the SWML estimates, the lines show the Schechter 
function fits. Note that one cannot fairly compare the shapes 
of the Schechter fits from this Figure, since the different 
redshift slices contain galaxies covering a different range of 
absolute magnitudes. In particular, the faint end slope of 
the highest redshift slice is very poorly constrained, since 
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Figure 3. Luminosity function for z < 0.1 (filled circles, con- 
tinuous line), 0.1 < z < 0.15 (open circles, dashed line), 0.15 < 
z < 0.2 (filled triangles, dot-dashed line), 0.2 < z < 0.3 (open 
triangles, dotted line). 
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Figure 4. Normalized comoving radial density plotted against 
redshift for the Northern (continuous line) and Southern (dashed 
line) Galactic hemispheres. Note that the radial density is nor- 
malized independently for the two hemispheres. 



there are no galaxies fainter than Afo.i r ~ —21.5 in this 
subsample. If analysis of the intermediate 0.1-0.15 redshift 
slice is limited to Mo.i r < —21.5, then the estimated faint- 
end slope changes from a = -1.17 ± 0.03 to -2.18 ± 0.44. 
The Schechter fits in Figure [3] are thus included for illus- 
trative purposes only, and we deliberately do not quote the 
Schechter parameters. 

The SWML data point for the lowest redshift slice at 
Mo.i r w —23 lies significantly above the Schechter fit. This 
is almost certainly due to errors in the photometry of the 
half-dozen or so objects in this bin. Of the ten apparently 
most luminous objects in the z < 0.1 slice (Mo.i r < — 22.6), 
all have the EDGE or C0SMIC_RAY flag set. Low redshift, lu- 
minous galaxies have much larger fluxes and apparent sizes 
than most galaxies, and hence are more susceptible to their 
measured fluxes being affected by cosmic ray hits or by lying 
near the edge of a CCD chip. 

Despite the uncertainties in the shape of the LF in 
the redshift slices, there is clear evolution in the amplitude 
of the LF, in the sense of an increasing amplitude (verti- 
cal shift) and/or luminosity (horizontal shift) with redshift. 
The estimated number densities of galaxies in the range 
—24 < Ma-i r < —21.5, obtained by applying @ to the 
non-parametric LF for each redshift slice, are given in Ta- 
ble For comparison, we also show the gala xy density in- 
ferred from the evolving parametric LF of iBlanton et alJ 
J2003d) . Their fit gives Schechter parameters a — —1.05, 
M * = -20.44 and 0* = 0.0149/i 3 Mpc" 3 at z = 0.1 in 
the 01 r frame, with characteristic magnitude M* evolving 
as M* = M * - Q(z - zo) with Q = 1.62. We ignore the 
negligibly small number density evolution in their fit and 
assume that the error in number density is dominated by 
the ~ 3% uncertainty in <f>* . Our density estimates agree 
within ~ 2(j for the three lower redshift slices, but the den- 
sity estimate for the highest redshift slice is ~ 9a larger 
than that inferred from the Blanto n et al. model. The re - 
sults presented here, and those of IBlanton et al] i2003d) . 
are entirely consistent, since the Blanton et al. parametric 
model assumed linear evolution of the characteristic magni- 
tude M* with redshift. Less than 6% of the galaxies in the 



full (—24 < Mo.i r < —16) sample are at z > 0.2, and so the 
Blanton et al. analysis would have been insensitive to rapid 
evolution at these redshifts. 

In the following section, we investigate evolution in the 
radial density of galaxies in narrower redshift bins. Note 
that the density estimator we use assumes a non-evolving 
luminosity function, so any effect found can equally well be 
ascribed to number or luminosity evolution. 



4 RADIAL DENSITY 

4.1 Estimator 

Just as the maximum-likelihood estimate of the luminos- 
ity function 4>(L) is independent of inhomogeneities in the 
galaxy distribution, one can also estimate the radial density 
p(z) independently of the assumed luminosity function by 
maximizing the likelihood 

C = T\p{zi) I [ maX( ' p(z)dV , (7) 

where z m i n (£i) and z m ax(£i) are the limiting redshifts at 
which a galaxy of luminosity Li would still be included in 
the survey. We fit p(z) by an arbi trary step function, u sing 
a variant of the SWML estimator dSaunders et aljll990l) . As 
in the maximum-likelihood estimate of 0(L), overall normal- 
ization is lost and so we have applied the constraint 

J p(z)S(z)w(z)dV I J S{z)w{z)dV = 1, (8) 

where S(z) is the selection function JSJ and w(z) is the 
weighting function @. This constraint is also used for the 
error estimates (cf. EEP). 

4.2 Results 

Our estimate of p(z) is plotted separately for the Northern 
and Southern Galactic hemispheres in Figure 2] We see a 
gentle increase in p(z) out to z ~ 0.2, with a steeper increase 
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Figure 5. Normalized comoving radial density plotted against 
redshift for all galaxies assuming the same cosmology and K- 
correction scheme as Figure [I] (continuous line) and assuming a 
A = cosmology (dashed line) and i<"-correcting to the r band 
(dot-dashed line). The open symbols show normalized densities 
in the four redshift slices used in the luminosity function analysis 
(Table □. 



to z ~ 0.3. This estimate of radial density is independent 
of the galaxy luminosity function only insofar as there is no 
correlation between redshift and luminosity. The observed 
increase in radial density could thus either reflect an increase 
in the number density of galaxies at higher redshift (whether 
due to number density evolution or the existence of a large 
local "hole"), and/or an increase in luminosity with redshift. 
If there is a large local underdensity, it is extremely unlikely 
that we happen to lie exactly at the centre of it, and so the 
consistent radial dependence of density in the Northern and 
Southern hemispheres provides strong evidence against the 
local hole hypothesis. 

We have also checked the dependence of the estimated 
radial density on the assumed cosmological model and our 
method of applying ^-corrections. Figure|^|shows estimated 
radial density for the full sample (north plus south) for 
the same assumed cosmology and if-correction as Figure [3] 
along with the estimate separately assuming a flat, A = 
cosmology, and A"-correcting the r-band galaxy magnitudes 
to redshift zero instead of redshift z — 0.1. The acceler- 
ated expansion in the A-dominated cosmology gives rise to 
a larger density at high redshift compared with the A = 
cosmology. Even though the increase in radial density with 
redshift is rather less if one assumes a A = cosmology, 
particularly for z > 0.2, the increase in density is still sig- 
nificant. The estimated radial density is entirely consistent 
whether one isT-corrects galaxy magnitudes to the 0,0 r or 
01 r bands, suggesting that any errors in the if-correction 
have a negligible effect on the estimated densities. 

We have tested the significance of the radial density 
evolution over the redshift ranges 0-0.15 and 0.15-0.3. In 
the lower redshift range, evolution is small but so are the 
error bars. In the higher redshift range, where evolution is 
more pronounced, the error bars are much larger. For each 
redshift range, we calculate the mean radial density and per- 
form a x 2 test of the null hypothesis that all bins have this 
mean density. We apply the same test to two non-evolving, 



Table 2. Results of a \ 2 test for uniform density within two 
redshift ranges (low-z: 0-0.15 and high-z: 0.15-0.3) for the DR1 
sample and for clustered and random simulated catalogues. In 
each case there are 14 degrees of freedom. 



Sample 


low-z 


high-z 


DR1 


594 


1730 


Clustered simulation 


400 


32 


Random simulation 


14 


11 



simulated catalogues. The first consists of a random distribu- 
tion of galaxies with compara ble galaxy numbers and L F to 
DR1. The second consists of a lSoneira fc Peeblesl l|l97Sf) hi- 
erarchical simulation with comparable galaxy numbers, LF, 
clustering and sky coverage to DR1. The resulting \ 2 statis- 
tics are given in Tabled in each case there are 14 degrees 
of freedom (15 data points minus one free parameter: the 
mean density). 

For the random simulation, the \ 2 values are consistent 
with the uniform density null hypothesis. 

For the clustered simulation, the \ 2 values are inconsistent 
with uniform density, particularly in the low-z range. This 
is to be expected given the presence of large scale structure 
in the simulated galaxy distribution. The biggest contribu- 
tion to \ 2 comes from the first three radial bins (z < 0.03) 
where density fluctuations are particularly pronounced since 
clustering is not smeared out by peculiar velocities in these 
simulations. 

The DR1 sample yields x 2 values marginally larger than 
the clustered simulation in the low-z range, and significantly 
larger than the clustered simulation in the high-z range. De- 
viation from a non-evolving p(z) is thus clearly very signif- 
icant in the redshift range 0.15-0.30 but only marginally 
significant in the redshift range 0-0.15. 

4.3 Comparison with LF 

The open symbols in Figure|S]show the densities of luminous 
galaxies inferred when normalizing the luminosity function 
in four redshift slices (Table 0, rescaled to have unit mean 
density. The horizontal bars attached to these symbols de- 
note the redshift range (the symbol is centred on the mean 
redshift in each bin, rather than at the centre of each bin) 
and the vertical bars denote the estimated error in density. 
The agreement between these two completely independent 
methods of estimating radial density is striking: the appar- 
ent increase in estimated radial density with redshift can 
be explained entirely by evolution of the galaxy luminosity 
function, with no need for a local underdensity. 



5 GALAXY NUMBER COUNTS 

In this section we explore the effect of galaxy evolution, as re- 
flected in our estimate of radial density, on counts of galaxies 
as a function of apparent magnitude. In par ticular, galaxy 
numb er counts in the APM Galaxy Survey jMaddox et alJ 
1990) are steeper than one would expect unless there is sig- 
nificant evolution of galaxies at low redshift. Can the evolu- 
tion seen here at z < 0.3 explain the steep APM counts? 
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Figure 6. Galaxy counts (per square degree, per unit magnitude) 
as a function of r-band magnitude. The continuous line shows pre- 
dicted counts assuming a non-evolving luminosity function, the 
dashed line assumes a radial density variation due to evolution as 
estimated in Section |3] Observed co unts are shown by di amonds 
(EDR), triangles (DR1) and circles iMaddox et al.lll990t APM). 
For the APM counts, we have applied the very rough correction 
from bj to r-band magnitudes r ~ bj — 1.3. 

The predicted number counts per unit magnitude are 
given by 

f°° dV 
n(m) = / 4>[L(m, z)]p(z)——dz, (9) 

Jo 

where L(m, z) is the luminosity of a galaxy at redshift z and 
with apparent magnitude m, p(z) is the radial density and 
dV is the comoving volume element at redshift z. These pre- 
dicted counts are plotted in Figure|S|for the non-evolving LF 
of Section f3 .3 . 1 1 with p{z) = 1 (continuous line) and for the 
same LF but with p(z) as estimated for the full sample in 
Section 2] (dashed line). We make the conservative assump- 
tion when evaluating © that density does not evolve further 
beyond a redshift of 0.3, so that p(z) = 2.0 for z > 0.3. 

For comparison, the open circles show number counts 
measured in the bj p assband from the APM Galaxy Survey 
iMaddox et al.ll 990*l , where we have made a very rough cor- 
rection to the r-passband using r « bj — 1.3. The evolving 
model derived from SDSS data is in remarkable agreement 
with the APM counts considering the different passbands 
used and areas of sky observed. If the steep APM counts are 
due to a local underdensity, then a very similar underden- 
sity exists in SDSS galaxy counts which come predominantly 
from the Northern sky. 

Observed galaxy number cou nts from the SDSS Early 
Data Release llYasuda et al.l l200Lj are shown as diamonds 
(we have summed the Yasuda et al. counts in the Northern 
and Southern equatorial stripes). These counts are signifi- 
cantly shallower than the APM counts and are consistent 
in shape with no galaxy evolution. However, these counts 
come from an area of only 440 square degrees and so are 
susceptible to fluctuations from large scale structure. 

We have therefore also estimated galaxy number counts 
from the 2099 square degrees of SDSS DR1. Galaxies are 
selected from the DR1 database according to the following 
criteria: 



(i) None of the SATURATED, BLENDED, BRIGHT or EDGE bits 
be set in the r-band photoFlags. 

(ii) The object be classified as a galaxy in at least two of 
the g, r and i bands. 

These counts are shown as triangles in Figure [(J While not 
quite as steep as the APM counts, they still lie closer to the 
evolving LF model than to the non-evolving model. 

Rapid evolution of the luminosity and/or density of the 
galaxy population at redshifts z < 0.3 thus provides a nat- 
ural explanation for the observed radial density of SDSS 
galaxies and steep number counts of APM and DR1 galax- 
ies. 



6 CONCLUSIONS 

We have presented evidence for significant evolution in the 
luminosity function of r-band selected galaxies in the SDSS 
DR1 at redshifts z < 0.3. This evolution gives rise to a factor 
~ 3 increase in inferred galaxy density betwee n redshifts 
and 0.3, in agreement with earlier findings from lEaiesI lll993l) 
and lBlanton et al] (12003d) . Folding this evolution into pre- 
dictions for galaxy number counts gives a remarkably good 
match to the slope of number counts observed in the APM 
Galaxy Survey and in DR1 itself. 

The aim of this short paper has simply been to demon- 
strate that the luminosity function of galaxies has evolved 
significantly at recent epochs, since z < 0.3, and that such 
evolution is sufficient to explain the steepness of the ob- 
served number-magnitude counts of APM galaxies, without 
the need to invoke a local underdensity or a magnitude scale 
error 1 . 

Note that these results do not preclude the existence 
of significant density fluc tuations in the loc al Universe on 
very large scales. Indeed, iFrith et al] <l2003l) have recently 
used galaxies in the 2 Micron All Sky Survey to indicate the 
presence of a region in the Southern Galactic hemisphere ~ 
200/i _1 Mpc in extent with a mean underdensity ~ 30%. This 
underdensity alone, however, is insufficient to explain the 
steepness of the APM counts. We have shown that evolution 
of the galaxy LF, possibly in combination with large density 
fluctuations, can explain the steep counts. 

Here we have not attempted to investigate the type 
of evolution that is occurring: the possibilities include any 
combination of luminosity evolution, density evolution and 
change in shape of the galaxy LF. Future work will use SDSS 
data to study galaxy evolution in detail. We plan to inves- 
tigate luminosity evolution in different passbands and for 
different galaxy types, and to perform a detailed investiga- 
tion of the spectral evolution of galaxies with redshift. 
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